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Abstract 

We show that nonlinearly elastic plates of thickness h —^ 0 with an e-periodic structure such 
that £~^h —>■ 0 exhibit non-standard behaviour in the asymptotic two-dimensional reduction 
from three-dimensional elasticity: in general, their effective stored-energy density is “discontin- 
uously anisotropic” in all directions. The proof relies on a new result concerning an additional 
isometric constraint that deformation fields must satisfy on the microscale. 


1 Introduction 

Understanding the behaviour of elastic plates (and, more generally, thin elastic structures) from 
the rigorous mathematical point of view has attracted much attention of applied analysts over the 
recent years. The related activity was initiated by the paper [5] by Friesecke, James and Muller, 
concerning plate deformations u with finite “bending energy”. This work appeared alongside the 
thesis [5D] by Pantz and was followed by a study of other energy scalings [S]. It puts forward the 
idea that homogeneous plates of thickness h, viewed as three-dimensional noirlinearly elastic bodies, 
afford a special compactness argument for sequences of deformation gradients Vu with elastic energy 
of order 0{h^) as h —^ 0. This argument is based on a new “rigidity estimate” [8] concerning the 
distance of local values of Vm from the group of rotations SO(3). The observation that this bounds 
the distance from a constant rotation field (which depends on u) allowed the authors of [8] to show 
that the “limit” elastic energy functional as h —>■ 0 is given by 

Eximiu) ■= <92(n(M)), U e (1.1) 

UJ 

where w C represents the mid-surface of the undeformed plate, II(m) = (V''«)^V'(9i'u A d 2 u), 
V' = { 81 , 62 ), is the matrix of the second fundamental form of an isometric surface m : w —>■ 
and Q 2 is a quadratic form derived via dimension reduction from a quadratic form appearing in the 
process of linearisation of elastic properties of the material in the small-strain regime. 

From the point of view of applications to real-world materials, it seems reasonable to ask in what 
way the above result is affected by a possible inhomogeneity of material properties of the plate in 
the directions tangential to its mid-surface. One can imagine, for instance, that the plate has a 
periodic structure of period e > 0 and try to replace the e, h-dependent family of h“^-scaled energy 
functionals by an “effective” functional in the sense of variational convergence [6]. Such analysis 
should reveal, in particular, whether the behaviour of the plate really depends on the relative orders 
of smallness of the parameters e and h. Work on this programme was started in the thesis m 
by Neukamm, where a series of similar questions was addressed for a periodically inhomogeneous 
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rod, i.e. a one-dimensional analogue of a plate. Further, a recent paper [14] has investigated the 
behaviour of periodic plates in the cases when h ^ e and h ^ e. This was followed by the work [22] . 
where the case e is addressed. In the present paper we develop an approach (see Sections 

H ED that in our view simplifies the derivation of the corresponding F-limit, via a “smoothing” 
approximation procedure that precedes the two-scale compactness argument, see Section Oof the 
present work. Smoothing is known to be useful in the asymptotic analysis of sequences of solutions to 
parameter-dependent PDF (see e.g. [10], [24]). In the present work we exploit similar considerations 
in the asymptotic analysis of bounded-energy sequences, where the energy is represented by an 
integral functional. Our smoothing approach replaces the approximation result m Lemma 3.8], 
when £ or h ^ e, and an additional technical statement (22] Lemma 3.7], when <C h <C e, by 
a two-scale compactness theorem for the second gradients of smooth approximations (Theorem 13.41 
below). Remarkably, in the situation where the deformation Uh € has no second derivatives, 

the smoothing procedure takes us to a setting where the second derivatives exist and are bounded 
in the L^— norm, see Lemma [221 

The added value of our approach, however, is revealed through its ability to deal with the more 
difficult case h <C which has remained open until now. An additional key observation in handling 
this case is that the determinant of the matrix of second derivatives of the smoothed approximation 
can be rewritten in a form amenable to the use of a compensated compactness argument (see proof 
of Theorem ih.l]) in order to derive an additional, ”fast-scale”, isometry constraint on the limit finite- 
energy deformations, see (na). It also perhaps worth mentioning that the smoothing approach yields 
a somewhat shorter route to the statement of compactness of finite-energy sequences in the case of 
homogeneous plates considered in [8]. 

Notably, a recent paper m contains an analysis of “zero-thickness” plates, where functionals of 
the “limit” form dia with explicit e-periodic a:-dependence of Q 2 are studied in the limit £ —>■ 0. 
Our results are consistent with those of nn, in the sense of offering an alternative route to the same 
“supercritically-thin plate” limit functional. However, while the authors of HU work with a two- 
dimensional formulation from the outset, our derivation is in the spirit of “dimension reduction” 
from the full three-dimensional problem. In addition, our approach offers a shorter route to the 
proofs of HH: which is checked directly. 

All of the above earlier works follow a traditional methodology of T-convergence relying solely on 
convergence properties of sequences with bounded energy combined with ’’structure theorems” for 
combinations of their gradients, without analysing their asymptotic structure as a two-scale series in 
powers of the parameters e, h. Our approach, on the contrary, is guided by an analysis of the original 
sequence of energy functionals via two-scale power series expansions. (For a detailed discussion of the 
viewpoint provided by the method of asymptotic expansions, see Appendix B.) For “supercritical” 
scalings <C it suggests in particular that, to the leading order, the approximating (“exact” or 
“recovery”) sequences have to satisfy an additional constraint 

det(II-h = 0, (1.2) 

where "0 = ip{x',y), x' € uj, y € Y := [0,1)^, is a term that appears in the two-scale limit of the 
original sequence of deformations and is responsible for the behaviour of the plate on e-scale. One 
key aspect that facilitates this approach is that the limit stored-energy function Q 2 is quadratic with 
respect to the second fundamental form of the surface, which makes the analysis of e, h-dependent 
plate energies somewhat amenable to a perturbation technique. 

The observation that the space of deformations with finite bending energy acquires the constraint 
(O, combined with a characterisation of isometric embeddings by Pakzad m and the earlier ansatz 
of [5], see (I7.76L allows us to carry out a rigorous proof of variational convergence. Thanks to the 
two-scale asymptotic structure, the presence of the second gradient of the corrector ip in the two- 
scale limit of deformation sequences provides a natural construction for the recovery sequences, cf. 
([4.3611 for the moderate regime. The supercritical regime /i <C is studied in Section [T] which 
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contains our main result, Theorem 17.21 the limit energy functional is given by 

Kira = ^/ Qrom(n(«)), 

UJ 

where 

Qhom = min J <32(2/, n + \ 7 y'tlj)dy, V' £ H‘^{Y) periodic, subject to (II. 2p . 
y 

The isometry constraint (11.21) implies, in particular, that is “discontinuously anisotropic in 

all directions of bending” as a function of the macroscopic deformation gradient Vu, see Theorem 
O and the subsequent discussion. To the best of our knowledge, this is a new phenomenon for 
nonlinearly elastic plates. Finally, we remark that the analysis of the “critical” scaling /i ~ is 
currently open. 

Throughout the text Xd denotes the d x d identity matrix, (2=2,3. Unless indicated otherwise, 
we denote by C a positive constant whose precise value is of no importance and may vary between 
different formulae. We use the notation di, i = 1,2,3, for the partial derivative with respect to Xi, 
dy^, i = 1, 2, for the partial derivative with respect to t/i, and Vy for the gradient with respect to y. 

2 Setting of the problem 

Let a; be a bounded convex^ domain in and let := w x (—f) be a reference configuration 
of an undeformed elastic plate, where h is a small positive parameter, 0 < <C 1. A deformation u 
of the plate is an iJ^-function from ilh into K.^. We assume that the material properties of the 
plate fl/j vary over w with a period e, which goes to zero simultaneously with h. For a plate made of 
a hyperelastic material (see e.g. [3]), the elastic energy of a deformation u is given by the functional 

I W(e~^x',Vu)dx, 
nn 

where x' := {xi,X 2 ) £ w. The stored-energy density function W{y,F) is assumed to be measurable 
and periodic with respect to y G Y := [0,1)^ and to satisfy the standard conditions of nonlinear 
elasticity (see e.g. my- 

A) The density W is frame-indifferent, i.e. W{y,RF) = W{F) for any R € SO(3) and any 

F G where SO(3) is a group of rotation matrices in 

B) The identity deformation u(x) = x is a natural state, i.e. W{y^X^) = mini?’ lU(y,F) = 0. 

C) There exists a constant C > 0 such that for aW y gY, F G he inequality 

lU(y,F) > Cdist2(F,SO(3)) (2.3) 

holds. 

Additionally we assume that (c/. [3], [5], [T3]) 

D) The density W, as a function of the deformation gradient, admits a quadratic expansion at 
the identity, i.e. there exists a non-negative quadratic form on such that 

lU(y,l 3 + G) = Q 3 ( 2 /,G) + r(y,G) (2.4) 

for all G G where 

r(y, G) = o(|G|^) as |G| —>■ 0 uniformly with respect to y GY. 

^The convexity of a; is a technical requirement which allows us to use the results describing the properties of 
isometric immersions for the supercritical case h e^. From the point of view of the aim of the present work this 
requirement is insignificant, since the limiting effective behaviour of an elastic plate is a local property. 
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We next rescale the transverse variable X 3 in order to work on a fixed domain. Namely, we 
multiply the transverse variable by h~^ so that the new variable belongs to the interval I := ^)- 

We reassign the notation X 3 to the new variable (this should not cause any misunderstanding since 
we will always deal only with the rescaled domain in the future), so a; € fl := a; x /. We also use 
the rescaled gradient 



f_d _^ 1 d \ 

\ dxi ’ 8 x 3 ’ h 8 x 3 ) ’ 


so that the rescaled functional is given by 


Eh{u) 



x', 'S/hu)dx, 


n 


where u G We consider sequences of deformations Uh in the bending regime, that is such 

that 

Eh{uh) < Ch\ (2.5) 

(In the non-rescaled setting this corresponds to the energies of order h^.) In this paper we focus on 
the case when h <C e. We consider e as a function of h, i.e. h is our “main” parameter and e is the 
“dependent” parameter. 

The non-degeneracy property (12.31) and the assumption (12.51) imply that 

J dist^ {VkUh, SO(3)) dx < Ch^. (2.6) 

n 

In what follows, we say that a sequence satisfying (1216]) has finite bending energy. 


3 Two-scale compactness and second gradients 

The two cases when e is of order h and h ^ e were considered in M- In these regimes the behaviour 
of the mid-surface of the plate is essentially macroscopic and the microscale only plays a role away 
from the mid-surface. In our case the microscopic behaviour of the mid-surface of the plate is very 
important. This is due to the fact that for each e-cell the corresponding piece of the plate itself 
behaves like a miniature plate. This “micro-plate” behaviour becomes even more dominant when 
h e^. In order to study this property we employ the method of two-scale convergence. The 
structure of the limit functional in the case of a homogeneous plate, namely, the fact that it is 
basically defined on the second gradient of the limit surface, suggests that the information about the 
microscopic behaviour should also appear in the form of a second gradient. The problem is that the 
original functional is defined only on iJ^-functions, and how one would obtain the second derivatives 
with respect to the fast variable y in the limit is not clear at the outset. However, while this is not 
possible in general, in the bending regime deformations possess an additional property which in a 
certain sense is equivalent to the existence of the second gradient, which we explain in detail in what 
follows. 

Let us recall some results from [8]. We denote 

Ya,h := a+ hY, a G hi?, uh := |J Ya,h- 

Yc,hCcj 

Theorem 3.1. Suppose that a sequence Uh from H^(rt) has finite bending energy. Then: 
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1. Up to a subsequence 


^hUh —>■ (V^u, n) strongly in 

where u = u(x') belongs to the space of isometric immersions, i.e. such maps u : uj ^ 

that = I 2 , and n := u^i A u, 2 ; 

2. There exists a piecewise constant map Rh ■ uJh ^ SO(3) such that 

j \VhUh-Rh\^ dx <Ch^] 

UlhXl 


3. The difference between the values of Rh in each pair of neighbouring h-cells is small in the 
sense that 

J max \Rh{x' + C) ~ Rh{x')\^dx' < Ch'^, 

i^h 

where consists of those ( G {—h, 0, h}"^ for which Fo+c,/i C LOh for the value a G hi? such 
that x' G Ya^h- 

Apart from the strong compactness of the gradients VhUh, which is true in our setting as well 
as for homogeneous plates, this theorem implies that oscillations of the gradients on the h-scale are 
bounded. The third property in Theorem 13.11 implies that the difference between values of Rh in 
neighbouring h-cells divided by h (which basically is the difference quotient) is bounded on average 
and a similar statement can be formulated for VhUh due to the second part of the theorem. Hence 
one can try to mollify Uh, expecting that the second gradient of the mollification will be bounded. 
In what follows we always assume that Uh is a sequence with finite bending energy. 

Let tp G C^{Bi), > 0, be radially symmetric, where Bi is a unit disc centred at the origin, 

p = 1 , and for each h > 0 denote by ph{x') := h~‘^p{h~^x'), a;' G the corresponding mollifier. 
Let uj' be a domain whose closure is contained in w. For each sufficiently small h > 0, consider the 
function Uh = Uh{x'), x' G tv', defined as the result of the simultaneous mollification of Uh with ph 
and averaging with respect to the variable xp. 


Uh{x') := j J ph{^-x')uh{f,X3)dfdx3= J j ph{i)uh{f, + x',X3) dfdx^, 


X G oj. 


(3.7) 


(The above integral is well defined in ut' for sufficiently small h since in this case the support of 
Phif — x') is contained in uj.) It is well known that the mollification Uh is infinitely smooth. We will 
need to estimate in the difference between Uh and the average of Uh with respect to X 3 , 

Uh{x') := J Uh{x)dx3, x' G uj\ 

I 

as well as its derivatives up to the second-order. Such estimates are obtained in the next statement. 

Theorem 3.2. Let w' be a domain whose closure is contained in uj. The mollification satisfies 
the following inequalities with a constant C > 0 independent of uj' and h: 


\\uh- Uh\\HJ{u:') <Ch, (3.8) 

W^'^UhWmu.') < c, 

< C. (3.9) 
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Proof. In the subsequent estimates we use a generalised Minkowski inequality (also sometimes re¬ 
ferred to as the Minkowski integral inequality): 


/ 


F{x,y)dnx{x) 


p 

dyy{y) 

i/p 

[/ 

F{x,y) 

V 

dyy{y) 


J 

X 



- 


i/p 


dfixix), 


where X and y are measure spaces with measures y-x (x) and y,y (y) respectively, F : d:” x —>■ K. is 

measurable and 1 < p < oo, see e.g. m- We obtain the L^-estimate in (13.81) as follows: 


\Uh — Uh\\L2{u') 


< 


{J {J +x') - Uh{x')) df^ dx'^ 

uj' 


R2 LJ' 


= max 
\i\<h 


f Vuh{ti + x')-£,dt < max \^\ f W'^'uhit^ + x')\\^^,,. 

Jo L^{u;') Jo ^ ’ 

< C'/i||V'u;i(x')||^2(^) < Ch. 


dt 


Now we estimate the difference of the gradients of Uh and Uh- First notice that by the properties of 
mollification, Vuh equals to the mollihcation of Vuh ■ Second, one has 

||V'm/i — Rh\\L'^{uj') < W^'uh — RhWl'^^q.') < Ch, (3.10) 


where Rt stands for the first two columns of the matrix Rh and Q,' := uj' x I. The bound (13.101) is a 
simple consequence of application of the generalised Minkowski inequality and the Holder inequality. 
Further, the following estimate holds: 


L^u') 


< 


\\Vdh-Vuh\\L-^^-) = I / g,h{£.){^'uh{^ + x')-Vuh{x'))d^\ 

B2 

/ Ph{i){Vuh{i + x')-V'uh{x')) df<max\\S/'uh{i + x')-V'uh{x')\\^^,,. 

J ^ ^ (o.ii) 

R2 

< max||v'u/j(^ -b x') - Rh{x' + Oh 2 /,, + max||.R/,(^ -b x') - Rh{x')\\ 2 /. 


-b||.R/i(a;') - V'u;i(a:')||^ 2 (^,) < Ch, 


where we use the generalised Minkowski inequality and Theorem 13.1 
Next we estimate the second gradient of Uh: given by 


= J - x')) 0 VuhiO d^ = - J ® ^'Mx' + 0 dO 

R2 B2 

where Viph{0 = h~^Vip{h~^0)■ Since the mollifier p is radially symmetric, d^^iphiO is an odd 
function of and an even function of ^ 2 , and analogously for hence we write 


d,V'uh = - J + -0)d0 * = 1,2. 

m\<h,ii>o} 
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Using the generalised Minkowski inequality once again we have, for i = 1,2 


{\i\<h,ii>0} 

< Ch~^ max||V'u;i(x' + 0 ~ V'm?i(x' — < C. 


L^u,') 




The last inequality is obtained analogously to (13.111) . 

Finally, we argue that V'u/j is bounded pointwise independently of h\ 

\V'uh{x')\= j iph{0'^'uh{^ + x')d^ < + 


R2 


< Ch W'V'ufii^ + x') — Rh{^ + C!h ^ C'l 

for all points x' G oj and sufficiently small values of h. 


(3.12) 


□ 


As one can see from the above, V'm?i is a “reasonably good” approximation of V'uh- It remains 
to approximate h~^dzUh by a smooth vector function. Let rih be the unit vector orthogonal to the 
surface Uh'- 

rih := diUh A d2Uh- 

We show that the L^-norm of Uh — Rh ,3 is at most of order h, where Rh,i stands for the i-th column 
of Rh- To this end we first transform the above expression for rih, as follows: 

n-h = diUfi A d2Uh = {Rh,l + {dlUh — Rh,l)) A {Rh,2 + id2Uh — Rha)) 

= Rh,3 + {diUh — Rh,l) A Rh,2 + Rh,l A {d2Uh — Rh,2) + {dlUfi — Rh,l) A {d2Uh — Rh,2)- 

Since Rh and Vuh are T°“-bounded, the L^-norms of the last three terms are of order h, hence 


\\nh-Rh,3\\L^iu.')<Ch. (3.13) 

Additionally, it is clear from (13.121) and (13.9p that 

l|V'n^|U2(^q < C. (3.14) 

Putting together the above estimates for rih, Theorem 13.21 and Theorem 13.11 we obtain the fol¬ 
lowing result. 

Proposition 3.3. Let uj' be a domain whose closure is contained in w. Under the assumptions of 
Theorem no the following estimates hold: 

11 (V Uh I Uh') V h'^h — Ch, 

11 (V Uh I ^/i) Rh 11 2^2 — Ch. 

As is customary in periodic homogenisation, and due to the fact that our problem is of two- 
scale nature, we require a suitable notion of convergence which takes into account the microscopic 
behaviour, e.g. the so-called two-scale convergence, introduced by Nguetseng in m- Out of several 
equivalent definitions of two-scale convergence, we will use the one based on periodic unfolding, see 
e.g. [5]. For convenience of the reader we give the definition of two-scale convergence that we use 
and some of its properties in Appendix Next we prove a general statement about the two-scale 
convergence of a second gradient, without any reference to a specific problem. This assertion may be 
used in various contexts where two-scale convergence of a second gradient is required. In the present 
paper we make use of this statement in Theorems 14.21 and 17.21 In the remainder of this section we 
consider an arbitrary d-dimensional space, returning to the case d = 2 in Section [5] onwards. 
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Theorem 3.4 (Second gradient and two-scale convergence). Let a; C be a bounded domain and 
Y := [0,1)'^. Assume that a sequence G H‘^{uj) is bounded in H'^{uj). Then up to a subsequence 

A v2^;(x)+V2^(x,y), (3.15) 

where v G H^{uj) is the weak limit in H^{uj) (hence, strong limit in H^(uj)) of the sequence and 
(fi G Lf {ui\ with zero mean with respect to y gY . 

Note that here Vg is a scalar function. The theorem directly applies to the vector case component¬ 
wise. 

Proof. Following the idea of the proof of the two-scale convergence for the first gradients in [5] we 
introduce the following function of x and y. 

Ve{x,y) := £~^{e~^{Te{Ve) - MviTeive))) - y^ ■ My {Te{VVg))^ , (3.16) 

where Tg is the unfolding operator (see Appendix El) . MY{f) := fy fdy is the operator of averaging 
with respect to the variable y, and y^ := y — (1/2,..., 1/2)^. We calculate the first and second 
gradients of Vg with respect to y. 

VyTg = £-\rg{yvg) - MyiTgiVvg))), Vpg = TgiV^Vg). 

Notice that the function Vg and its gradient VyVg have zero mean value in y. Hence, we can apply 
the Poincare inequality: 

— C'll^y'*^e||L^(i.iJxV) — ^*11 ^ C, 

where the last bound follows from the assumptions of the theorem. Hence, the sequence Vg is 
bounded in L‘^{oj-,H'^{Y)') and it converges weakly in this space, up to a subsequence: 

Vg ^ tip weakly in (w; iJ^(T)). 

Let a function tp G Lf{uj-, H^{Y)) be such that 

^ ^ y"" • (V^u)y°-b yj. (3.17) 


Clearly, under this notation one has 

NyVg V^u(x) + Vytp{x, y) weakly in Lfipj x Y), 

which is equivalent to (13.151) . Now it only remains to prove the periodicity of tp with respect to y. 
To this end, let ip{x, y') be an arbitrary function from C^(uj x Y'), where y' := (yi,..., y^-i)^ and 
Y' := [0,1)'^“^. We write the difference between values of Vg on the opposite faces of the cube Y 
corresponding to y^ = 1 and ya = 0 (tbe proof for other components of y being exactly the same), 

e{x,y',l) -Vg{x,y',0) = £~'^ {Tg{vg){x,y',1) - Tgivg){x,y',0)) - MyiTgiddVg)) = 

= £~‘^{Te{vg){x,y', 1) - Tg{vg){x,y',0) - Myidy^Tgivg))^ = 

= (%{ve){x, y', 1) - Tg{vg){x, y, ^) - J {%(yg){x, y', 1) - Tg{vg){x, y', 0))dy'). 

y 
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Notice that Te{ve){x, y', 1) = Te{ve){x + ecdi 0), where ea is a unit vector in the direction of Xd 
axis, hence we derive that 


j {ve{x,y',l)-Ve{x,y',0))iljdxdy'= 


ojxY' 


£ ^(%ive){x,y',0) - j %{vs){x,y',0)dy'^{tlj{x - £ed,y') - ilj{x,y')) dxdy' 


(3.18) 


ujxY' 


The sequence £ ^{'tp{x — eed,y') — tp{x,y')) converges to—9^'!/'strongly in x T'). Let us denote 

Ze{x,y) ■■= e~'^{Te{ve){x,y) - J Te{ve){x,y)dy'). 

Using a slight variation of the argument in Proposition 3.4 in [5] one can see that converges 
strongly in L'^{uj;H^{Y)) to (y“)' • V'v as £ —>■ 0. Hence, its trace 

Ze{x, y', 0) = £~^{Te{Ve){x, y',0)- J Te{Ve){x, y', 0)dy') 

converges to (y°)' • Vv weakly in x Y'), see Proposition 13.51 below. Passing to the limit in 
(13.181) yields 

lim J {ve{x,y',l)-ve{x,y',0))iljdxdy' = - j {y^)' ■ V'v ddilj dxdy'. 

uixY' uixY' 

We carry out the the same calculations for z{x, y) := ^ y'^ • (V^u) y°, which results in 

z{x, y', 1) - z{x, y', 0) = (y'')' ■ ddVv. 

Multiplying the last expression by ■;/' and integrating by parts we get 

j {z{x,y',l) - z{x,y',0))tpdxdy'= - J {y’^)'■ V'v ddi’dxdy'. 

ujxY' ujxY' 

Comparing the latter with (13.191) we conclude that y:>{x, y) is periodic with respect to y. 
Analogously we prove the periodicity of 'S/yip. Indeed, we have 

VyUe(a;,y',l) - VyVs{x,y',0) = £“^ (7;(Vue)(x, y', 1) - 7;(Vue)(a;, y', 0)). 

Note that by the assumptions of the theorem and properties of two-scale convergence (see Appendix 
the sequence %{Vve) converges to Vv strongly in L^{uj-, H^{Y)). Hence, by Proposition 13.51 
below the trace Te{Vve){x,y',0) converges to Vv weakly in L'^{uj x Y'), and therefore 

lim J (yyV^{x,y',1) - VyV^{x,y',0))tpdxdy' 

ujxY' 

= lim J %iyve){x,y',0) £~^{'il){x - eed,y') -'il^{x,y')) dxdy'= - J Vv dd^) dxdy'. 

ujxY' ujxY' 

Since Vyz{x, y', 1) — Vyz{x, y', 0) = ddY7v we have 

J {Vyz{x,y',l)-VyZ{x,y',0))tpdxdy'= - J Vv ddi^ dxdy', 


uixY' 


ujxY' 


which implies the periodicity oiVyp. Summarising all of the above we conclude that p G (w; iJpgj. (U)). 

□ 
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Proposition 3.5. Suppose that a sequence fn{x, y) converges to f{x, y) weakly in H^{Y)) and 
strongly in LF‘(uj x Y). Then its trace on uj x dY converges to the trace of f weakly in Lf(ui x dY). 

Proof. Without loss of generality we will assume that / = 0 and prove the assertion for the part of 
the boundary corresponding to yn = 0. We need to show that 

lim 

n—¥oo 

ujxY' 

for all ijj G C^{uj X Y'). Let us fix -0 and consider it as a function defined on w x P via 
if{x,y) := ilj{x,y'). Let us define function hn{y) as J^fnf’dx. It is easy to see that /i„ € H^{Y) 
and \\hn\\L2(Y) < CWfnWmuxY), \\hn\\m{Y) < C'll/nIIi(Y))• Hence hn ^ 0 weakly in H^{Y) 
and strongly in Lf{Y). By the trace theorem —>■ 0 strongly in L^{Y'). In particular, 

Jy' dn{y', 0) dy' 0. But 0) = /„(a;, y', O)'0(x, y') dx, hence (13.201) follows. □ 

Now, we have all necessary tools to proceed to the L-convergence statement. 


fn{x, y', O)0(a;, y') dxdy' = 0, 


(3.20) 


4 Limit functional and F-convergence for the “moderate” 
thickness-to-period ratio e 

Our main result is establishing L-convergence of the rescaled functional Eh to the limit homogenised 
functional defined on isometric surfaces. L-convergence is a standard tool to describe the asymp¬ 
totic behaviour of nonlinear functionals, introduced by De Giorgi [B]. The general definition of 
L-convergence is the following; 

Definition 4.1. Let Fi, i = 1,2,... and iGo be some functionals defined on functional spaces TLi 
and TLoo- Then we say that the sequence of functionals Fi V-converges to the functional F^, 

Foo = L - lim Fi, 

i—¥oo 

if the following two properties hold: 

1. Lower bound. For any sequence Ui G TLi converging in a certain sense to Uoo G TLoo one has 

liminf Fi(iti) > Foo{uoo). 


2. Recovery of the lower bound. For any Uoo G TLoo with Foo{uoc) < oo there exists a 
sequence Ui G TLi converging to Uoo in a certain sense such that 

lim Fi{ufj Fx:(^oo)- 
2^00 

In the above definition the notion of convergence depends on the underlying properties of the 
functionals and may be quite non-standard, in particular, because the spaces TLi and TLoo often 
do not coincide, as in the present paper, for example. Moreover, the notion of convergence in the 
second part of the definition may be different from the convergence in the first part. In this case 
the former must normally be “stronger” that the latter in order to facilitate the main purpose of 
the definition (nevertheless, this requirement is not always critical, e.g. see [2]). The lower bound 
inequality ensures that the limit functional gives in the limit a lower bound for the values of the 
functionals T). The recovery property implies that the lower bound Foo is the greatest lower bound. 
In particular, the definition implies that 

lim (inf Fi) = inf Foo, 

2—>00 
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and, if a sequence Ui of “almost minimisers” of Fi converges to Moo, then the latter is a minimiser 
of Foo- 

In the present work we obtain different F-limits for the regimes e ^ h ^ and ^ h. In 
this section we present our results for the “moderate” regime e ^ ^ e^, for which the limit 
homogenised functional is given by 

Korain) := ^ / Ohom(n)dx', M € HUcu) , 
u 

where II is the matrix of the second fundamental form of u, 

II := n = diu A 92M, 

and the homogenised stored-energy function is given by 

<9hom(II) := min / Q 2 (y, H + (4 21 ) 

Y 

where <52 is the quadratic form defined on 2 x 2 matrices and is obtained from the form Qa by 
minimising when reducing to two dimensions (notice that Qa vanishes on skew-symmetric matrices): 

Q 2 [y, A) ■= min Q 3 (?/, ^ -I- 63 (g) a -I- a (g) 63 ) = min Q 3 {y,A + a0 63 ). (4.22) 

aGR^ aGM^ 

where we implicitly assume that the space matrices is naturally embedded in the space 

by adding zero third column and zero third row, cf. [ 8 ]. 

Theorem 4.2. Suppose that e ^ h ^ e^, i.e. e~^h = o(I) and h~^e^ = o(l) as ft. —>■ 0. Then 
the rescaled sequence of functionals h~‘^Eh T-converges to the limit functional in the following 

sense. 

1. Lower bound. For every bounded bending energy sequence Uh S ift^(O) such that 'S/'uh 

converges to V'm weakly in u € F{?^^{uj), the lower semicontinuity type inequality holds: 

limmi h~^Ehiuh) > £'i™^(u). 

h^O 

2 . Recovery of the lower bound. For every u € there exists a sequence 

such that converges to V'u strongly in Lfljl) and 

\iTah-‘^Ehiu]f^) = £^“om(M). 

R —>-0 

Proof. Lower bound. Since WhU is close to a rotation-valued function Rh we can use the frame 
indifference property of the stored energy function W and its Taylor expansion near identity to 
linearise the functional. Roughly speaking, we will use the observation that 

W(y,Vf,u^) = W(y,R/^Vf,Uf,) = W(y,I 3 + (R/lV- 13 )) « Q 3 (y,hGh), (4.23) 


where 

G,, := ft-i(i?^V^M,,(x)-l3)- 

However, the last asymptotic equality in (I4.23P has to be made precise, since the smallness of the 
term hGh is not pointwise. We also need to understand the structure of the (weak) two-scale limit 


II 




of Gh- To this end we will use a special representation of Uh- Let uj' be a domain whose closure is 
contained in u> and fl' = uj' x I. For small enough h we can write Uh in fl' in the form 

Uh{x) = Uh{x') + hx^nhix') + hzhix), 

where Uh is the mollified mid-surface of Uh defined in the previous section and rih is the corresponding 
normal. Then one has 

Gh = h-^Rj{(y'uh\nh) - Rh) + xsRjiVuhlO) + RjVhZh- (4.24) 

Notice that Gh is bounded in by Theorem 13.11 The first term on the right-hand side is 

bounded due to Proposition 13.31 The second term on the right-hand side is bounded by (13.141) . 
Hence, the third term is also bounded. 

Since we perform the dimension reduction when passing to the limit, we are interested only in 
the behaviour of the top left 2x2 block of the matrix Gh- Let us use the following notation: if 
A G then A denotes the first two columns of A, and A is obtained from A by omitting the 

third column and the third row. We have 

Gh = h-^Ri{Vuh- Rh) + X3RlV'nh + RjVzh- (4.25) 

The limit of the first and the third terms on the right-hand side of (14.251) is of no importance for 
us apart from the fact that they do not depend on X 3 . This is obvious in case of the first term. As 
for the third term, let AlyxiiTeizh)) ■= consider 

Wh ■■= £~^{Te{zh) - AlYxl{Te{zh))). 

It is easy to see that S/yWh = Te{Vzh) and dsWh = d^e~^Te{zh) = e~^hTe(h~^d‘iZh)- Hence we have 
||Vy,x 3 W/i||LP(n'xV) < C. Since Myxi^Wh) = 0 we can use the Poincare inequality ||w/i||LP(n'xV) < 
G\\S/y^x 3 Wh\\Lp(n'xY)- So Wh is bounded in LP{uj'\W^’'p{Y x I)). Then, up to a subsequence, we 
have 

Wh w{x,y) weakly in LP{^uj';W^’P(Y x /)). 

On the other hand we know that d^Wh vanishes in the limit, hence w is independent of X 3 : 

w = w{x',y). (4.26) 

Therefore, we have 

h-^Rj{Vuh-Rh)+RlVzh^Hix\y) (4.27) 

for some H G L^(w' x Y), up to a subsequence. 

Now we examine the expression {Rh)^Vnh- The element (1,1) of this matrix can be written as 


Rh,i ■ diUh = Rh,i ■ diiUh A d2Uh + Rh,i ■ diUh A di2Uh = 

= diiUh ■ d2Uh A Rh,i + di2Uh ■ Rh,i A diUh = 

= —diiUh ■ Rh,3 + diiUh ■ {d2Uh — Rh,2) A Rh,i -I- di2Uh ■ Rh,i A {diUh — Rh,i), 

where the vector product operations take priority over the inner products. Consider the second term 
in the last expression (the third term can be dealt in the same way). By Theorem 13.II and Theorem 
O one has 


\\dllUh\\L^(uj') < C'l \\d2Uh — Rh,2\\L^{u:') < Ch' \\d2Uh — 7?^,21| < G. 

On the one hand, the second term is bounded in L^(w'), and so converges weakly two-scale to some 
function from x Y) up to a subsequence; on the other hand, for any G^(uj' x Y) test-function 
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we get zero in the limit. We conclude that the weak two-scale limits of the second and the third 
terms are zero. We can apply analogous argument to the other elements of {Rh)^Vnh. Therefore 
we see that {Rh)^ynh and —\7'^Uh ■ Rh ,3 have the same weak two-scale limit (the “dot” product 
gives a 2 X 2 matrix whose elements are —dijUh ■ Rh, 3 )- Applying Theorem 13.41 to Uh one has 

V'^Uh ^ +Vlipix\y) 

up to a subsequence, for some (/? G {VL'\ , whereas Rh ,3 converges strongly to n due to 

Theorem 13.II Hence, 


{RhVV'nu ^ ( - V'^u + Vl^p) -71 = 11 + Vl'iP, 


(4.28) 


where tp := —+ ■ n 
subsequence 


G L^(fl';i?pgr(F)). Putting together (I4.27|) and ()4.28p we conclude that up to a 
Gh ^ X 3 (ll(x') + V2 V>(x',2/)) + H{x',y). (4.29) 


We can proceed to the last step of the proof. In order to use the Taylor expansion of W we need 
to control Gh pointwise. Let Xh be a characteristic function of the subset of on which hGh is 
relatively small, namely. 


Xh ■■= 


1 ii\hGh\<h^/^, 
0 otherwise. 


(4.30) 


It is easy to see that /q,(I — Xh)dx —^ 0. Since W is non-negative we can restrict the integral to 
and discard the set of points on which \hGh\ > h}^'^ by using Xh'- 


= h 


-2 


h ‘^Eh{uh) = h J W{e ^x',yhUh)dx > h J W{e ^x',XhyhUh)dx 

n n' 

e~^x',Xhil3 + hGh))dx = j Q 3 {e~^x',XhGh)dx + h~'^ J r{e~^x',XhhGh)dx. 


(4.31) 


By the assumptions on W we have r(e ^x',XhhGh)\ < giXh\hGhf) for some function g{t) = o{t). 


An elementary lemma below implies that 

J r{e~^x', XhhGh)dx ^ 0. 


(4.32) 


Lemma 4.3. Suppose that t~^g{t) —>-0 as t —> 0, and let {/(} C L^{Xl') be a family of non-negative 
functions such that t~^ ft is bounded with respect to t and |l/i||Loo(n') —>• 0 as < —>■ 0. Then the 
convergence 

J 9 {ft{x))dx ^ 0 


holds. 

Proof. Notice that 

J 9{ft{x)) 


dx = 


giftjx)) ftjx) 
ft{x) t 


dx < 


sup 


9 {s)\ f Mx) 


^s6(0.||/dL=(Siq) S 


dx —0. 


Here it is implied, with some abuse of notation, that {ft{x)) ^g{ft{x)) = 0 whenever ft{x) =0. □ 
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We apply Lemma 1431 to the expression (14.3211 by setting t := h?, ft := Xh\hGh\^- Notice that 
the particular choice of the cut-off threshold in (I4.30p is not important: instead of one can take 
any threshold that goes to zero slower than h. 

Now we estimate the quadratic term in (14.311) . By definition of the form Q 2 ant the properties 
of the unfolding operator we have 


J Q'i{e ^x' ,XhGh)dx> J Q^ie ^x',XhGh)dx > j J Q2{y,%{xhGh))dydx ( 433 ) 


O' Y 


It is well known that convex integral functionals are lower semi-continuous with respect to the weak 
convergence in functional spaces, see e.g. [7]. Assume that we deal with a weakly converging (as in 
(|4.29l) ') subsequence Gh- Then passing to the limit and recalling that Xh converges strongly to 1 we 
obtain the following inequality: 

lii^inf J J Q 2 {y,TeixhGh))dydx > J J Q2{y,X3{ll{x') + Vl'ip{x',y)) + H{x',y))dydx 

O' Y 

= ^jj Q2(2/,n +-b J J Q2{y,H)dydx (4 34 ) 

uj' Y uj' Y 

-hj J Q2iy^^ + ^l^)dydx' >^l I Q^omimydx'. 


O' Y 


Note that in the second step in the above we integrate with respect 0:3 € I, which causes the cross 
term of X 3 (II -f V^'0) and H in the quadratic form Q 2 to disappear. At the last step we just used 
the definition of 

A simple argument by contradiction shows that the inequality 


lii^inf J J Q 2 {y,Te{XhGh))dydx > ^ J J QZmi^i)dydx' 


(4.35) 


O' Y 


ui' Y 


holds for the whole sequence Gh- Putting together (I4.31L (I4.32|l . (14.3311 and (14.3511 and noting that 
w' is an arbitrary subset of ui concludes the proof of the lower bound. 

Remark 4.4. The above proof applies to the general case of h e, up to the last inequality in 
/I4-341)- IPe will use this observation in our analysis of the lower bound in the case h <C in Section 

0 


Recovery sequence. In order to build a recovery sequence one normally needs to ensure that 
the argument of the stored-energy density W is an function, so that one can pass to the limit 
using the Taylor expansion (12.41) . If this condition is not in place one would have to estimate W 
on the “bad” set where the argument might be uncontrollably large, which would require additional 
restrictions on W. For example in ([ 8 ]) the authors used a truncation result for Sobolev maps, 
which consists of changing a function on some subset to make it together with its gradient essentially 
bounded and estimating the measure of this set. In this paper we will use a result of Pakzad from [12] 
(which appeared after [ 8 ]). It states that the space of infinitely smooth isometric immersions from 
a convex set ZU C into is strongly dense in the Sobolev space of isometric immersion 
Let u G (~l C'°°(aJ) and II be the matrix of the corresponding second fundamental form. Let 

also tjj{x', y) and d{x'^ y) be correspondingly a scalar and a vector functions from G°°{oJx Y) periodic 
with respect to y. 
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Denote := ipi^x'^x'), := d{x',e ^x') and consider the sequence 

2,2 

Uh = u — e^ipen + hxsTie + h^-^Rds, (4.36) 

where n := diuAd 2 U, rie := n + edy-^tpediu + edy 2 'ip^d 2 U and R := (V'^jn). Calculating the rescaled 
gradient of Uh we obtain 

'^hUh = (V'u - eVy{'ipen)\n^) + hxsiV'n + Vy^dy^^tpediu + dy^4!ed2u)\Rde) + Uh, 
where Uh denotes the rest of the terms: 

Uh -.= -e'^{V{Pj,n)\0) + ehx3{V'{^y,^Pe^lU + ^y,^|;,^2u)\0) + + e-^yy){Rd,)\0). 

Due to our assumptions we see that Uh = O(e^) = o{h) uniformly on uj. A simple calculation shows 
that 

(V'm - e\7y{'tpen)\ns)^{V'u - e\7y{'ipen)\ne) 

( dy^tpedy^lpe 0 \ 

=13 + e^ dy^lpedy^lpe {dy^'^ef 0 \=l3 + o{h). 

\ 0 0 J 

(4.37) 

This implies that there exists an infinitely smooth function Rh with values in SO(3) such that 

Rh {V'u - eVy{ipen)\ne) = la + o{h). 

(It can be constructed explicitly via Gram-Schmidt orthonormalisation for the columns of the matrix 
(V'u - eVy{ipen)\ne).) 

Further, consider the expression 

Gh := h~^{RlVhUh -T3) = X3RJ{V'n + Vy{dy^i!ediu + dy^4ied2u)\Rde) +o(l). 

Notice that Rh converges to R uniformly on w, converges strongly two-scale to and dg 

converges strongly two-scale to d. Hence, the convergence 

Gh ^xsill+ VlPj\d) (4.38) 

holds. Passing to the limit in 


h '^Eh{uh) = h ^ 


W{£ ^x',VhUh)dx 


^x',Gh)dx + h ^ 


r(e ^x',hGh)dx 


a 


n 


n 


we see that the remainder term vanishes since Gh is uniformly bounded in L°°(uj), and for the 
quadratic term due to (14.381) we have 

J Q 3 {e~'^x',Gh)dx = J Q 3 {y,Te{Gh))dxdy + J Q 3 {e~^x',Gh)dx 

n ClxY AeX/ 

^ ^ / Qa {y, (n + Vli;\ d)) dx'dy 

uixY 


(for the definition of Ag see Appendix |A|. Thus 


lim h '^Eh{uh) = ^ [ Q 3 {y,{ll + '^l'ip\d))dx'dy. 

h-i-O 12 I a / 

UIXY 


(4.39) 
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We conclude the proof with a standard argument showing that the above assumption of smooth¬ 
ness of M, tp, d is not restrictive. Let u be an arbitrary element of and take 'ip{x',y) € 

Hp^^{Y)) to be a solution to the minimisation problem in (14.211) corresponding to the matrix 
of the second fundamental form of u, that is 

Qhon.(n) = J Q2{y,ii + vlpj)dy. 

Y 


Let also d{x\y) be such that 


Q 2 ( 2 /, II + vg= Qa (y, (II + Vl^\d)). 


There exist sequences of C°° functions uj, tpj and dj converging to it, ip and d in the corresponding 
function spaces such that 


J Qsiv, (IIj + Yyipj I dj))dx'dy - J QZn,i^i)dx' 

ojxY uj 



For each j define via formula ()4.36p with u, 'ijj and d replaced by Uj^ iljj and dj. For each 
sequence the analogue of convergence (14.3911 holds. It follows that there exists a subsequence 
hj —>• 0 such that 

Ehiul^^)- I QZnmdx' 

U) 

Hence it™'^ := is the required recovery sequence. □ 



5 A microscale isometry constraint on the limits of finite 
bending energy sequences in the “supercritical” regime 

h 

During the process of constructing a recovery sequence we have encountered a problem in the case 
when h <C e^. Namely, introducing the corrector terms —e'^ip^n and edy^ipediu + edy 2 iped 2 U in (14.361) 
incurs an error of order (c./. (|4.37l) ') which becomes of order h~^e^ when scaled in the subsequent 
argument. This error does not vanish if h or h ^ e^. In order to understand the nature of 

the problem we write a formal asymptotic expansion for a recovery sequence, see Appendix [B] In 
the process of eliminating the errors of order and higher, it transpires that ipe needs to satisfy 
some solvability condition which is equivalent to the zero-determinant condition (15.401) below. As we 
show in the next theorem, this condition is intrinsic for the limits of finite bending energy sequences 
in the case h <C e^. This is due to the fact that II -|- approximates the matrix of the second 
fundamental form of the mollified mid-surface Uh of the plate (defined in Section [3|) on e-scale. It 
appears that Uh has to be almost isometric in the case h <C e^ for the expression {V'uhZV'uh —I 2 
to be of order o{h). Noting that the matrix of the second fundamental form of an isometric surface 
has zero determinant suggests the condition det(II-I-= 0, where ip is the two-scale limit of the 
sequence ipe- 

Theorem 5.1 (Zero determinant condition for the two-scale limit isometric surface.). Let Uh he 
a finite bending energy seguence, be the mollified mid-surface of the plate defined by and 
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oj' g bj. Assume that h <C . Then the weak two-scale limit of in L‘^(ui' x Y), 

u € € L^{u!'; Hp^^(Y)), satisfies the “zero determinant” condition 

det(II + Vlf)) = (IIii + dl^y^il)){ll 22 + dy^y^fj) - (IIi 2 + = 0, (5.40) 

where tf '■= ~T ' 

Proof. The proof is inspired by the general idea of the proof of the identity det II = 0 for the 
isometric surfaces (see [H]). In our setting, however, the mollified mid-surface Uh of the plate is 
not exactly isometric. This, together with the two-scale nature of our problem, requires in-depth 
analysis on the microscale. 

The following identity is easily checked by a straightforward calculation (recall that Uh is a C°° 
function): 

deiV^Uh = dl2{diUh ■ d2Uh) - ^{d22idiUhf + d‘i^{d2Uhf)- 

Note that we can consider the third order tensor as a 2 x 2 matrix whose elements dijUh are 

vector-valued functions. Thus when we take its determinant we understand the multiplication of 
entries as the scalar product of vectors. We will use this convention henceforth. We can rewrite the 
above identity applying the unfolding operator: 

TeidetV'^Uh) = (dy^y^Te{diUh ■ d2Uh) - ^{dl^TeidiUh)'^ d-dy2Te{d2Uh)'^)y (5.41) 

We consider the left-hand side of the last identity as a scalar product of two “vectors” fe := 
{Teid^iUh),Teidf 2 Uh)V and ge := (Teid^ 2 Uh),-Teld^ 2 Uh)V■ Since curly/e = 0 and divy = 0, 
we can apply the compensated compactness theorem, see e.g. |23j . which implies that the left-hand 
side of (15.411) *-weak converges to the scalar product of their weak limits, that is 

lim / / Te(detW'^Uh)r]dydx' = / / detiV'^u d-Vlip)r]dydx' 

^^0 J^, Jy Jui' Jy 

for any g G (u>'xV). Note that x' plays the role of a parameter, and does not affect our argument. 
Let us now consider the right-hand side of (15.411) . Due to Proposition l3.3l we have diUh = Rfi,id-0(h), 
h —>■ 0 , i = 1 , 2, where 0 (h) is understood in terms of the L^-norm. In particular, diUh - d 2 Uh = 0{h), 
[diUhY = 1 + 0{h). Thus 

■ d2Uh) - ]^{dl(re{diUhf P dl2re{d2Uhf)^ 

= ^~\dly20{h) - i(a^yi + 0{h)) + dl2{\ + 0(/i)))) (5-42) 

= 9ly20{hle^) + dl2p{h/e^) + 52.0(V£^). 

Multiplying the last expression in (I5.42|) by an arbitrary function g G C^{bj' x Y) and integrating 
by parts twice we see that it converges *-weakly to zero given that h/e^ = o(l), i.e. 


f f det{y'^u{x') d- Vyip{x', y))gdydx' = 0 . 

Juj' Jy 


This implies 

det(V'2M(a:') -f ^pix', y)) = 0. (5.43) 

At the last step of our proof we need one additional property of the two-scale limit (p{x' ,y). We 
know that the second-order derivatives of u{x') are parallel to the vector n{x') : one has 


dijU ■ Rk = diRj ■ Rk = djRi ■ Rk = di{Rj ■ Rk) - Rj ■ dkRi 


-Rj-dkRi=d, i,j, fc€{l, 2 }. 
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since at least two of the indices i,j, k coincide. Notice further that 


—dijU ■ n 


dijU ■ {d 2 U A diu) + d^jU ■ [diU A i9im) = diU ■ {diju A d 2 u) + diU ■ [diu A d 2 ju) 
= diU ■ dj{diu A d 2 u) = diU ■ djU = 11^-, i,j € {1, 2}, 


(5.44) 


hence dijU = —Ilyn. In particular, one has det = — det II. In order to derive (15.4011 from ()5.43p 
we need a similar property for ip to be held. In fact, we require a much weaker property of dy^ip 
being orthogonal to Ri for i = 1,2. We establish it with the help of the following statement. 

Proposition 5.2. Let 

Ze{x,y) := {Te{uh) - MyiTeiuh))) - MviTsiV'uh)) 

(cf. iS.16]) ). Then 

\\dyiZe ■ 'Te{Rh,i)\\L^{ui' XY) 0 OS h —?► 0, i = I, 2. (5.45) 

In the proof of this proposition we will need the following simple inequality, which is a sort of 
Poincare inequality. 

Lemma 5.3. Let f G H^{Y). Then 


I 1/(2/) - fiOl^dydC < 2 1\Vf{y)\^dy. 

YxY Y 


Proof. Since the set of C°°{Y) functions is dense in H^{Y) we only need to prove the lemma for 
/ G C°°{Y). First we write down a pointwise estimate for \f{y) — /(C)|- 


1/(2/)-/(C)l 


\/f-dl 


'Cl 


j‘V2 

dif{t,C, 2 )dtP / 82 f {Cl, t)dt 

JC2 


< 


< [ (|v/(t,C2)| + |v/(Ci,t)|)dt. 

Jo 


We finish the proof with 


J \f{y)-f{C)\^dydC< J j\\Wf{t,C 2 )\^ + \^f{Ci,t)ndtdydC = 2 j\Wf{y)\^dy. 

YxY YxY Y 


□ 


Proposition \5.2[ Taking the gradient of the function 2 ^ and using the properties of the unfolding 
operator we have 

VyZ, = e-\Te{Vuh) - MviTeiVun))). 

Let us write the L^-norm of the right-hand side (dropping the coefficient e~^) in a special form: 


\\%{Vuh) - MF(re(V'n„))||i.(^,^^) = J TeiV'uh) - J Te{V'uh)dC 

ui'xY Y 


dydx = 


w'xV 


{Te{V'uh){x, y) - Te{V'uh){x, C))dC 


Y 


dydx < 


(5.46) 


< J \Te{^'uh){x,y) -Te{^'uh){x,C)f dCdydx 

ui'xYxY 
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Applying Lemma [531 to 7^(V'u/i) we have 


cj'xYxY 


\%{V'uh)ix,y) -%(y'uh){x,Of dCdydx < 2 \\Vy%(Yuh)\\l 2 ^^,^Y) = 

= 2e^%iV'^dn)\\h^^,^Y) < C'e^ 


(5.47) 


where in the last inequality we used the property of boundedness of the second gradient of Uh, see 
Theorem l3.2l Due to Proposition 13.31 we can replace Te(Vuh) in the last estimate by Te{Rh)' 


J Te{Rh){x,y)-Te{Rh){x,0 dCdydx <Ce^. 

cj'xYxY 

Similarly to (15.461) we can write the expression in (15.451) as 


(5.48) 


J e ^{Te{diUh){x,y)-Te{diUh){xX)) ■Te{Rh,i){x,y)dCdydx. 

uj'xYxY 

As follows from Proposition 13.31 in order to prove the statement it is enough to show that 
J {Te{Rh,i){x,y)-%{Rh,i){x,C)) ■%{Rh,i){x,y)dCdydx ^ 0, i = 1,2. 

cj'xYxY 

Notice that the integrand in the latter is non-negative everywhere. The expression under the sign 
of integral has the form (a — 6) • a at every point, where a and b are unit vectors. We use the 
following elementary identity, 2(a — b) ■ a = |a — h\^, together with (I5.48P to conclude the proof of 
the proposition: 


J {Te{Rh,i){x,y)-Te{Rh,i){x,0) ■Te{Rh,i){x,y)dCdydx 

uj'xYxY 

= J \Te{Rh,z){x,y)-Te{Rh,z){x,0\'^dCdydx <Ce. 

uj'xYxY 


□ 

Theorem l3.4l implies that the sequence Zg converges weakly in (w'; H‘^(Y)) to the function (cf. 

iHD) 

^ + yty2did2U + ^{y^fd^u + ip. 

Since Te{Rh) converges to R strongly in L^{uj' x Y), we have 


dy,Ze ■ Te{Rh,i) ^ dy,z- Ri, 

*-weakly, i.e. with test functions taken from C^{uj' x Y). (In fact, with some effort it is possible to 
show weak convergence in x Y), however, *-weak convergence is sufficient for our purposes.) 

Further, since the second-order derivatives of u are parallel to the vector b, we have 


dyi Z. * Ri — ip * Ri. 

Lastly, it follows from Proposition 15.21 that the *-weak limit of dy-Ze ■ Te{Rh,i) equals zero, hence 

dy^ip • i?i = 0, 1 = 1,2. 
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The above orthogonality property implies that the derivatives of ip can be represented in the form 
dyi^P = {dy^ip- R 2 )R 2 + {dy^ip-n)n, dy^ip = [dy^tp- Ri)Ri + {dy^ip■ n)n. Differentiating further we get 
the following identities for the second derivatives of ip (in particular we see that the mixed derivative 
is orthogonal to both Rx and i? 2 ): 


-^ 2)^2 + ■ n)n, 

= i ^ V 2 ^ ■ -^1)^1 + 

dy^dy^ip = {dy^dy^ip ■ n)n. 

Substituting the latter and V'^it = —Iln into (15.431) and recalling that if) = —ip ■ n we obtain 

(lOOl) . □ 

6 The structure of isometric immersions 

The idea of the construction of a recovery sequence for the supercritical regime h <C £^, which is 
performed in the next section, relies on the properties of isometric surfaces from There has 

recently been a number of works devoted to this topic, see m,m, [T5] , [19] . The results of Section 
|5|(as discussed at the beginning of Section |7|) imply a discontinuous dependence of the homogenised 
energy density function on the direction of the bending. This observation drastically reduces the 
number of manipulations one can perform with a surface in order to construct a recovery sequence. 
In particular, we can not use the result on approximation of surfaces by infinitely smooth 

ones because this approximation does not preserve the orientation of the bending directions. In 
our paper we will use mostly the results of m, which we review below. We then give an example 
of an isometric surface for which the set of rational directions is irregular (c/. Theorem 17.11) . and 
prove a statement on regularisation of isometric surfaces required for the construction of a recovery 
sequence. 

6.1 Auxiliary observations 

For any u € domain tu can be partitioned into “bodies” and “arms”, according to the 

following definition. 

Definition 6.1. A body is a maximal subdomain of ui on which u is affine and whose boundary 
contains at least three segments inside oj. 

A differentiable curve F : [0, Z] —>■ a; is referred to as a leading curve ifu is affine when restricted 
to the lines normal to F([0,Z]), which we refer to as leading segments. In the approximation 
procedure below we also use the term “leading curve” for the curve 7 := (u o F) on the surface. 

For a given leading curve F, generatrices are the lines tangent to the surface u at the points of 

7 and orthogonal to the latter. 

An arm u;[r] is a subdomain of lo covered by leading segments corresponding to a leading curve 
F. (See Remark \6.3\ below.) 

Remark 6.2. Note that a maximal subdomain of uj on which u is affine and hose boundary contains 
at most two segments inside oj is always an arm or part of an arm. 

Remark 6.3. Suppose that T is a leading eurve parametrised by its arclength, and denote by T := F' 
the tangent vector and by N := the unit vector orthogonal to F. Then the corresponding 

arm is given by 

a;[F] = {r{t) + sN(t) : s G e [0,;]} Dw. (6.49) 

For each t G [0, 1], the set {F(t) + sN{t) : s G M} D w is the corresponding leading segment. 
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Since leading curves are not uniquely defined, it is more reasonable to define the vector field N 
on the whole a;[r] as a continuous vector field of unit vectors parallel to leading segments, and the 
vector field T as the continuous unit vector held orthogonal to N. In particular, we have T = T' 
for any leading curve F. In what follows we consider N and T either as vector helds N{x'), T{x') 
or the normal N{t) and the tangent T{t) of the curve T{t), depending on the context. We will say 
that the value of T (equivalently, iV) is rational if the ratio of its components is rational. We will 
also say in this situation that the corresponding leading segment is rational. 

Remark 6.4. It can be shown that: 

a) For any given body or arm the part of its boundary that is situated in the interior of tv is a 
union of straight segments whose vertices belong to dux, 

b) As follows from all leading curves are in the class A combination of this observation 
with the result of m establishing that C (w) implies the fact that leading curves 7 on the 
surface u are in the class . 

Since the leading segments do not intersect in w, the vector held N (and, hence, T) is locally 
Lipschitz in w. In is not difficult to prove, using a geometric argument, the following lemma. 

Lemma 6.5. The local Lipschitz constant 


^ \Nix') - N{x'o)\ ,. \Tix')-Tix’o)\ 

Cl{xo) := hmsup , , _ ,, = limsup ——^ 


Xr, G UJ 


a;' — a;; 


X' — Xr. 


satisfies the estimate 

Cl{x'q) < dist“^(a:Q, duo). 

This lemma implies, in particular, local Lipschitz continuous differentiability of leading curves. 

Note that for any leading curve F, unit tangent vectors of 7 are also tangent to the surface 
u : —>■ K.^. Therefore, the surface u is reconstructed on each arm from the knowledge of the 

corresponding leading curve and the non-vanishing principal curvature. Indeed, denote by Kr the 
curvature of F, so that F" = T' = kyN- For a given 7 , we dehne a Darboux frame (r, v, n) by the 
tangent r = 7 ' to the curve 7 , the unit vector v = {V'u)N tangent to u and orthogonal to r, and 
n = T A i/. This Darboux frame satishes the system of equations (see e.g. m P- 277], [T9]) 


( 0 Kr ^ \ ( X 

—kt 0 0 I I 

—K 0 0 J \ n 

which has a unique solution subject to appropriate initial values. Given 'y(t) and u{t) we can write 
uL[r] as 

u = 7 (t) -I- sv{t), (6.51) 

where s parametrises the generatrix corresponding to the value t on the curve F. 

Remark 6.6. There is freedom in the choice of the sign of Kr and the direction of the vector field 
N. We make this choice in such a way that the normal n in the above Darboux frame coincides with 
the normal n defined as diu A d 2 U in the discussion preceding the eurrent section. 

Let us write down the formulae for for the gradients V'u, and the matrix II in terms of the 
leading directions and curvatures. Since x' is related to {t, s) via x'{t, s) = F(t) + sN{t). We have 




d{xi,X2) 

d{t,s) 


{T' + sN\N) 


Ti(I - sKr) -T 2 \ 

Tjfil-SKr) J 
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and therefore 


d{t,s) _ ( Ti(l-sKr) T^^{l-sKr) 
d{xi,X2) V “^2 


The above expression is well defined since 1 — SKr is positive everywhere on u;[r], cf. Lemma |6.51 
Then we can write the gradient and the second gradient of u in terms of the geometry of the curves 
7 and L: 


V'u = (Tit — T 2 V, T 2 T + Tiv ), 


v u = 


1 — SK-p 


n T 1 T 2 n \ 
T 1 T 2 n T^n J ' 


(6.52) 


Notice that {cf. (15.441) 1 


II = - 


K 

I — SKp 


Tf T,T 2 \ 

T 1 T 2 Tf ) ■ 


(6.53) 


Let us introduce the following notation 


( Tf T1T2 
V T1T2 r| 


(6.54) 


As will be shown in the next section, the homogenised stored energy density Q^om supercritical 

regime is in general a discontinuous function of its argument 11. More precisely, it is continuous with 
respect to the coefficient —k(I — sKr)”^, but discontinuous with respect to II (the direction T). 
This suggests that one can modify the curvature k in a continuous way {i.e. without changing much 
the homogenised energy), however, a perturbation of the vector field T may lead to uncontrollable 
changes in the homogenised energy. 


6.2 Regularisation of irregular surfaces. 

The results of [19] and m assert some regularity of isometric immersions u € Hf^^{uj), namely, 
u G C^{uj) and its associated leading curves in w have locally Lipschitz derivative. Nevertheless, 
isometric immersions can be quite irregular for the purpose of this paper. To illustrate this, we give 
an example of a surface for which the set of “rational” directions for the corresponding leading curve 
is the union of a fat Cantor set and a set of zero measure. 

Example. Let us recall the construction of a fat Cantor set. At the first step we remove from 
the middle of [0,1] an open interval Bi of length 1/4 and denote the remaining set Ai. At the next 
step from the middle of each the two subintervals of Ai we remove open intervals of length 1/16 
each, the remaining set is denoted by A 2 and its complement by B 2 and so on: at the n-th step 
we remove 2"“^ disjoint intervals of length 1/2^", whose union is a set of measure 1/2"+^. Thus we 
have Ai(A) = Xi{B) = i, where A := DnAn, B := By Xd we denote the Lebesgue measure 

in The Fat Cantor Set A is closed and is nowhere dense {i.e. contains no interior points). 
For each n we define a leading curve via its unit tangent vector T„(t) given as follows. Let ip be 
a function that is continuously differentiable on the interval [0,1], has zero mean and is such that 
V?(0) = if{l) = (^'(0) = p'{l) = 0, e.g. p{t) = (sin(27rt))^, and set Aq = [0,1], Bq = 0, Tq := (1,0)^. 
For each n = 1, 2,..., we define T„ := T„_i on the set [0,1] \ (i?„ \ B^-i), and on each component 
interval In^k of the set \ B^-i we define 

"Tra '■= ~ ^ k = 1,2, ..., 

where ipn,k = — tn,k) is the scaled (by a coefficient /3„, |/3„| < 1, and the length of the 

interval 1/2^") and appropriately translated (by the left end tn,k of the interval In,k) version of p, 
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which “fits” into this interval. We then integrate to obtain the curves 

t 

r„(t) := y r„, t e [o,i]. 

0 

Note that at each step of this procedure we “redefine” the vector r„ only on the set Bn \ Bn-i- 

We take /3„ = /3/2^”, where 0 < /? < 4(niax is a constant that will be specified later. So far 

we have a sequence of curves r„ with a bounded second derivativdl: F" < C/3. It has a subsequence 
that converges strongly in H^{0, 1) to a curve F G H^{0, 1) with |F"| < C/3. In particular, one has 
F' = (1,0)^ on the fat Cantor set defined above, and the subset of B on which F' takes rational 
directions has zero measure. Let w := [0,1] x [0,1], where / < 1 is defined by the right end of F, 
namely it is the first coordinate of the point F(l). If we choose /3 such that |F"| < C/3 < (diam(u;))“^, 
the lines normal to F do not intersect within w. (Clearly, one has diam(a;) < so the choice 
/3 < l/(C-\/2) satisfies the required condition.) Lastly, we can construct an isometric surface u via 
procedure described above (see formulae ()6.50|1 and (I6.51|l 'l with some “reasonable” curvature k, e.g. 
K=l. For thus constructed u the domain ui is the only arm and F is its leading curve in the sense 
of Definition 16.11 

The above example illustrates that for a leading curve T{t) the set on which T{t) has rational 
directions can not in general be represented as a union of intervals (up to a null set). In what 
follows we argue that one can replace an arbitrary isometric surface m by a “sufficiently regular” 
isometric surface, in order to obtain an explicit construction of a recovery sequence. Note that one 
can not simply approximate u by an infinitely smooth isometric surface, since general, 

is an everywhere discontinuous function of its argument. Thus the matrix II, see (I6.54|l . of the 
replacement surface should be different from the one of u only on a relatively small set. First we 
need to recall some properties of isometric surfaces. We reformulate and complement some results 
of [19] (Lemmas 3.8 and 3.9) in the following statement. 

Theorem 6.7. Let u G For any ^ > 0 there exists v G such that 

X 2 i{x' : u{x') ^ v{x')}) + ||u - v\\h2(u,) < <5, 

the surface v has a finite number of bodies and arms, each maximal connected region with non¬ 
zero mean curvature lies between two regions where v is affine, the corresponding leading curves are 
Lipschitz differentiable and do not have common points with dui. 

The above properties of v are partly derived from the proofs of Lemmas 3.8 and 3.9 in m- 
Next, we briefly elucidate them. Every maximal connected region (we denote it by E) where u has 
non-zero mean curvature is attached to either one or two regions where u is affine. If E is situated 
between two regions where u is affine, then the lengths of the leading segments that cover E are 
bounded below by a positive constant. Lemma 16.51 implies that on any closed subset of w the vector 
field T is Lipschitz continuous. In particular, we can cover E by a finite number of leading curves so 
that they do not have common points with the boundary of w, and hence have Lipschitz continuous 
derivative. If E has just one “affine” region attached to it or E = w, then it may happen that the 
infimum of the length of leading segments covering E is zero. Then one can not control the Lipschitz 
constant on the sequence of segments which lengths tend to zero (such sequence will converge to a 
point on the boundary of w). In this case we can replace m in a small region near the boundary point 

direct calculation shows that 

ir"! < 

which e.g. for the choice g>(t) = (sin(27rt))^, t G [0,1], and (3 < 2\/2 gives the value C = 4-71-^2/3. 
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where such degeneration of leading segments occurs by an affine function making sure that the new 
surface is still an element of Then we get to the situation discussed previously. 

As was mentioned at the beginning of this section we aim at the construction of recovery sequences 
for the homogenised functional in the supercritical case h <C given by (17.661) . Theorem l6.7l 
implies, in particular, that for ^ > 0 there exists v G with a finite number of bodies and 

arms that depends on S so that 

< S. 

The next theorem ensures that we can replace u G with an isometric surface such that 

the set on which its leading segments are rational can be represented as a union of finite number of 
arms and a set of arbitrary small measure. 

Theorem 6.8. Let u G E[?^^{uj) be such that the subset ujr{u) of u! covered by rational leading 
segments has positive measure. Then for each (5 > 0 there exists a surface us G such that: 

a) The surface us has a finite number of bodies and arms, whieh we generically denote by a;[r]; 

b) The subset ujr{us) of co covered by rational leading segments can be represented as 

uJr{u6) = USs, 

Fr 

where each arm w[rj.] consists of parallel rational segments (in particular, each leading curve Tr is 
a straight segment and M(a;[rr]) is cylindrical); 

c) The set Ss is measurable; 

d) The estimate 


M{Ss) + \\u-us\\m;^) + \E(0„,{u) - E(0^^{us)\ < 6 


(6.55) 


holds. 

Moreover, us can be chosen in such a way that on each arm the corresponding curvature k (cf. 
fd.JOl) and i6.53\} ) is infinitely smooth and 11s G where the latter in the matrix of the second 

fundamental form of us. 

Remark 6.9. The conditions that the curvature k is smooth and II 5 G L°°(w) is reguired for the 
construction of a recovery seguence. In particular, we take advantage of these facts when using the 
Taylor expansion formula for W, see i2.4\ ). 

Proof. In view of Theorem 16.71 we may assume without loss of generality that u has a hnite number of 
bodies and arms from the outset. Let us consider an arm a;[r] with the corresponding leading curve 
r(t), t G [0,/], and assume that the subset K = K{T) of [0,1] on which the directions of T = T{t) 
are rational has positive measure Ai(Ar) > 0. Since the set of rational directions is countable, we 
can order its elements in such a way that Ai(Ari) > Ai(Ar 2 ) > ■ • where Ki is the preimage of 
a particular rational direction Ti, i = 1,2,3,..., so that K = UiKi. We take an arbitrary rj > 0 
and let m = m{r]) be such that for := \JfL^Ki we have Ai(Ar \ Krf) < p. Note that all Ki are 
closed, in particular, the distance between any two of these sets is positive. This implies that we 
can remove a finite number of open intervals Li C [0,1] \ i = 1... ,k, so that the remaining set 
:= [0, /] \ (J^^^Li) (we may also remove the points 0 and I from if they are isolated) possesses 
the following properties: 

1. The estimate Ai(T^ \ AT^) < 77 holds. 

2. The set is the union of a finite number of mutually disjoint closed intervals Lj, j = 1,... ,n. 

3. For each j one has Lj fl AT^ = Lj fl Ki for some i G {1, 2,..., k} and both ends of the interval 
Lj belong to Ki. 
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Let us consider a particular interval Lj = [oj, bj]. The part of w[r] covered by the restriction of 
r to Lj, which we denote by is contained between two parallel segments r(aj) + sNi and 

r( 6 j) + sNi corresponding to the rational direction T^. In fact, loosely speaking, u;[r|i^] is mostly 
covered by the leading segments which are parallel to Ni. On the rest of the set, corresponding to 
the restriction , where the leading segments are not parallel to Ni, we replace them with new 

leading segments that are parallel to Ni. Notice that the difference between N{t) and Ni on the set 
Lj \ Ki is not greater than Crj, where C is the Lipschitz constant of N{t). We can write 

w[r|Lj] = {r(t) + sNi, s G K.t G Lj} n w. 

We carry out this procedure for each interval Lj. In this way we change the leading segments 
covering a;[r] only on a small set We estimate the measure of this set as follows: 

A 2 (a;,,[r]) < Xi{Lri\ Kri)diain.{uj) < 77diam(w). (6.56) 

Let us denote by N^ and the vector fields corresponding to the new set of leading segments. We 
have 

iV,=iV,T, = Tinc.[r]\w,[r], 

ll-^r; — -^|U“(i.d,[r]) —>■ 0, \\Trj — 'r||L~(c^,,[r]) 0 aS Tj ^ 0. 

We denote by the union of all rational leading segments that do not belong to the set (J" =i^[rUJ 
and estimate its measure as follows: 

^ 2 {S^) < Xi{K \ A"^)diam(a;) < 77 diam(a;). (6.58) 

For each arm we build a family of leading curves F^ and isometric surfaces Ujj, which will be 
shown to converge to the surface u uniformly as 77 —?► 0. We assemble the approximating parts Urj 
on all arms and bodies in a continuous fashion. We then choose a sufficiently small value of 77 so 
that the resulting surface denoted by us satisfies the inequality ()6.55[) . 

We define F^ as the maximal solution in aj[r] to the following ODE: 

r;(t) = r,(r,(t)), F,(o) = r(o), (6.59) 

where T,, is the field of unit vectors orthogonal to the modified set of leading segments of a;[r]. Notice 
that F^ uniformly tends to F as 77 0. In particular, for small enough values of rj, the new curve will 

cover the whole a;[F], i.e. one has w[r] = w[r,,]. We denote the domain of F,, by [0, noting that 
in general 7 ^ /. Both T an are Lipschitz continuous on any compact subset of w, the measure of 
the set w,;[r] on which T 7 ^ is proportional to 77 . This implies that on w[F^] one has a point-wise 
estimate \T{t) — r^(t)| < Cr], where C is proportional to the Lipschitz constant. It follows that the 
solution of (16.591) converges to F uniformly as 77 —0. To each leading segment of w[F] correspond 
unique t G [0, 1] and i G [0, /,,] such that F(t) and r^(t) belong to this segment. Let us denote by Cr? 
the corresponding map : 1 1 —>■ t. It is not difficult to see that is differentiable with converging 
to I as 77 —0. A direct calculation yields the following expression for the directional derivative of 
T{x') in the direction of T{x'): 

where {t, s) is such that x' = F(t) -|- sN{t). In particular, on the set a;[F] \ a;^[F] we have 

oC,(t) ^ »r(t) as., ^ 0. (6.60) 

Now we can construct a new surface on w[r,,]. To this end we will use the equations analogous 
to (|6.50l) and (16.511) . First, we need to adjust the non-zero principal curvature of the surface k to 
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compensate the shift of the leading curve Fj, along the leading segments, so let 


K{t) 

We define the leading curve on the surface by 

i 

Iv^i) ■= J 7^7(0) = 7(0), 

0 

with a unique Darboux frame satisfying the following system of equations: 

/ T-'(t) \ / 0 ^ Kr^{i) KrjO(-^i) \ / Tr,{i) \ 

\ n'r,{i) J \-Kr,oC-^{i) 0 0 / \ nr,{i) J 

with the initial conditions (t^( 0), n^(0)) = (t( 0), i^(0), n(0)) , where is the inverse of Cry- 

Then the new surface Urj is defined by 

s) — T 

for suitable values of t and s. Substituting t = Cri{t) we see that (r^ o 1/^ o (r]{t),'rirf o is 

the solution to the system 

^ / T-ryOCryW \ f ^ O Kr^it) \ / Tr, O (^{t) \ 

— =C;^(t) _Kp^oC^(t) 0 ® Vr,oC,r,{t) , (6.63) 

\ n^oCvit) j \ -Kr;(0 0 0 J \ nr,oC^{t) J 

with the initial conditions (r^ o Ct7(0), i^ry o Cr;(0)) ° Ci;(0)) = ('^(O)r *^(0)> ^(0)) • Due to (16.6011 and 

the properties of Cr; the coefficients on the right hand side of (j6.63ll converge to the coefficients of 
system ()6.50p in L°° norm. This implies 

\\{T^o(^,iyr, O (ri,nn O C,,) - -)> 0, 

hryOCr; -7 ||l“[ 0.;]-t 0, (6.64) 

W^V ~ “||L~(<.,;[r]) 


Analogously we have 


V Lirj — {^r),l'^r) '^r),2^r)i 2'^7] H“ : 


7'^U, = 


1 - s«:r„ 


TliUn 


DjypTjy 2 

J- TJ.l J- rj 


t 77,1-^ 77,2 ^77 


^ry,2 ^ly 


where the terms on the right hand sides are functions of t G [0,4], cf. (j6.62p . Let us substitute 
i = Ciy(t) and s = s — (T,, o Cry(t) — r(t)) • N{t) in the above, so that r,;(£) + sNn{i) = r(t) + sN{t) 
everywhere in a;[r] \ a;,,[r]. We obtain, in particular, via (16.6011 and (I6.61|l that 


CHt) 


i{t) 


1 - sKrJt) 


= ■^ 1 - sKr(<) 

K{t) 


on a;[r] \ u;^[r], 
on w,, [r]. 
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Then from (I6.56p . (16.5711 and (16.6411 we derive the following convergence properties: 


\v Ur^ ^ ^llL“(cj[r]) 0 

V Urf-V '“IL2(^[r]) 0' 


V lAri — V W Tost rri\ rrn —^ 
II^ = II on a;[r] \ a;,,[r], as ry —>• 0, 


(6.65) 


where II,, is the matrix of the second fundamental form of Ut^ . 

The estimates (16.5811 . (I6.65P imply the validity of (16.5511 . Indeed, since the number of bodies and 
arms of u is finite, we can carry out the above construction for every arm (on bodies the surface 
remains unchanged), and then assemble an iJ^(a;) isometric surface starting from an arbitrary arm 
or body and attaching sequentially to a current piece the next one by applying a rigid motion, so that 
the resulting surface satisfies the properties analogous to (16.6511 . Then choosing a sufficiently small 
rj and denoting the corresponding surface by us and the union of all by Ss yields the estimate 
(j6.55p . Finally, we split each arm of us into smaller arms w[rr] corresponding to the cylindrical 
parts with rational directions according to the construction and their complement. 

We conclude the proof with an argument for additional regularity of us, which is claimed in the 
statement. On each arm we can approximate the corresponding curvature k in the L^-norm by an 
infinitely smooth function. We then reconstruct the pieces of the new surface using the equations 
(j6.50ll . (16.511) and 7 ' = t with appropriate initial conditions, and attach them together in a fashion 
similar to the above. Taking sufficiently good approximations of k we ensure the validity of the 
estimate (|6.55p for the new surface. 

The following scaling argument allows us to obtain the boundedness of II^ in w. Due to the 
smoothness of k and ()6.53l) the matrix of the second fundamental form II^ is bounded on any compact 
subset of oj. Without loss of generality we can assume that w contains the origin. Since w is convex 
by the assumption, it is a star-shaped domain, i.e. uirj := (1 — r])u! C ui, where we use rj again to 
denote a small positive parameter. Let us define the scaled surface: us^r]{x') := {l — ri)~^us(^{l — ri)x'^, 
x' G oj. It is easy to see that us^n is an isometric surface (this is why we need to use the coefficient 
(1 — r])~^ in front of us). In addition, one has Il 5 ,,,(a;') = (1 — r 7 )Il 5 ((l — ri)x') and II^,,, € L°°{uj). 
Therefore, for any quadratic form Q, in particular for obtain 

J Q{lls,^ix'))dx' = (1 - 7)2 J Q(ll5((l - r,)x'))dx' = (1 - 7 ) J Q{llsix'))dx'. 

(jJ (jJ LJjj 


Hence, choosing small enough 7 and re-denoting ,, by us we arrive at (j6.55p . 


□ 


7 Limit functional and F-convergence in the “supercritical” 
regime <C 


As follows from the compactness result of Section and the isometry constraint of Section [5l in the 
case when h <C we should define the homogenised energy functional as 


KoM :=^l Qron,(n)dx', 

(jj 


(7.66) 


with the elastic stored-energy function 


Qrom(n) := min J Q 2 (y, II + det(ll += 0 . (7.67) 
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Zero determinant condition is rather restrictive in the case of periodic functions. Next we state a 
theorem describing the structure of periodic functions if) subject to this condition, first recalling that 
II can be represented in the form (16.5311 . The coefficient —k(I — sKr)~^ is one of the two principal 
curvatures of m, the other principal curvature is always zero. 

Theorem 7.1. Let II be a non-zero matrix of the form i6.5S\} and let tp € be such that 

det(ll +=0. (7.68) 

Then the derivative of ip in the direction {—T 2 , Ti)^ is zero, i.e. one has 

fp{yi,y2) = tpriTiyi + T2y2) 

for some ipT- In particular we have ipT = V’t(?/i) or ipx = '0 t(2 /2) ond ipx is 1-periodic for T = 
(1,0)^ and T = (0,1)^ respectively. For other rational T such that T 1 /T 2 = p/q € Q, where p/q is 
an irreducible fraction, the function ipT is periodic with period P = |Ti/p| = \T 2 /q\. For irrational 
vectors T the function ipT (and, hence the function ip) is constant. 

Proof. The proof is based on a result by Pakzad in [H] . Let us denote 

'o{y) ■= + 2/12/2II12 + ^ylih 2 + Ip, 

so that det(VyU) = 0 is equivalent to equation (j7.68() . Then Vu is continuous and for any regular 
convex domain D in and any y £ D there exist a segment of line passing through y and connecting 
two points on the boundary of D such that VyU = lly + Vyip is constant along this segment. Since 
the periodicity and continuity (by the Sobolev Embedding Theorem) of ip imply its boundedness, it 
follows that the segment along which VyV is constant has to be parallel to the vector {—T 2 , Ti)^: 
along this direction Ily is constant. Hence ^yip is constant in the direction {—T 2 , Ti)^, and the 
derivative of ip in this direction has to be zero due to the boundedness of ip. The rest of the statement 
of the theorem easily follows from the periodicity of ' 0 - CH 

The above theorem implies that at the points x' where the principal directions of the isometric 
surface u are irrational we have ( 3 hom(II) = /y 02 ( 2 /, 11 )^ 2 /, while on the rational directions the 
corrector ip is non-trivial and in general one has (5 |qjjj(II) < Jy Q 2 {y,LV)dy. This means that in 
general Qhom(I^ everywhere discontinuous function of the direction T. We can write ip in the 
form 

0(221:2/2) = -:;— -0t(7"i2/i + 122 / 2 ), (7.69) 

1 — sKr 

where ipxit) is the solution of the minimisation problem 

min J Q 2 {y,Tlil+ip!f{Tm+T 2 y 2 )))dy, 0 T(t) £ (7-70) 

y 

with II defined by ()6.54p . Notice that the Sobolev embedding theorem implies that iprit) is contin¬ 
uously differentiable. 

Theorem 7.2 (T-convergence). Suppose that h , i.e. he~^ = o(l) as h ^ 0. Then the rescaled 
sequence of functionals h~^Eh T-converges to the limit functional in the sense that 

1. Lower bound. For every bounded bending energy sequence Uh £ such that 'S/'uh 

converges to W'u strongly in u £ the following inequality holds: 

limmf h~^Ehiuh) > E^^^{u). 
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2. Recovery of the lower bound. For every u € there exists a sequence G 

such that converges to V'm strongly in L^(fl) and 

R—>-0 


Proof. Lower bound. The proof follows the derivation of the lower bound in the case when 
e ^ h » up to the last inequality in (14.341) . whose right-hand side can have replaced by 

Qhom given by (17.671) . thanks to Theorem 15.II 

Recovery sequence. A natural method of approximating the deformation of an elastic plate 
is to first write an expression u = u(x') for the mid-surface of the deformed plate, which takes into 
account the microscopic oscillations in the definition of the homogenised functional (represented by 
the function f}), and then “add thickness” to this surface. One possible approach to deriving the 
mid-surface expression is via an asymptotic expansion, with u{x') as the leading order term followed 
by a series of correctors in sequential powers of e. In the case h <C e^, the first corrector involving the 
function ip is shown to satisfy equation (17.681) . see Section|Bl which is the solvability condition for the 
equation at the next order in e. However, higher-order equations all involve two unknown functions 
and their solution is not straightforward. As was already mentioned in the previous section, in our 
construction we adopt a different approach based on the properties of isometric immersion, which 
we describe next. The details of the formal asymptotic argument are given in Section [BJ 

As was shown in Theorem 1 7. 11 the corrector ip is constant whenever T is irrational. This implies 
that if the measure of the set on which T is rational is zero, then the limit elastic energy is simply 

^ y Qhom(n)da;'= ^ J Q 2 {y,U)dydx'. 

OJ LJXY 


Thus, the only case 
measure. 

Let u G H?^. By 
that 


that remains to be considered is when the above mentioned set 


Theorem l6.8l we can choose a sequence of regularised surfaces usj 

1 


L;ro„>(w)-i?ron.K) < 

' J 

J J Q2{y,lhi)dx'dy <^. 


Ss. Y 


has non-zero 
0, such 


(7.71) 


We first construct an approximating sequence of plate deformations for us^ for each j and then apply 
a diagonalisation procedure to obtain a recovery sequence for u. Notice that all objects below depend 
on the parameter 6j, however we drop it from the notation in most of the cases. Let us consider 
an arm wjTr] of us^ corresponding to a rational direction T. On this arm the leading curve L^ is a 
segment of a straight line, T is constant, Kr^ = 0, the principal curvature is —K{t) (see (16.531) 1. and 
the vector v of the Darboux frame (r, v, n) of us^ o T^ is constant. Let ip = —n{f)ipT{Tiyi -I- 222 / 2 ) 
be the minimiser of (I7.67|) . c/. also (17.691) . We have 


TiTa 

T1T2 T| 

Qrom(II)=min J Q2{y,ll-KVlipT{,Tiyi+T2y2))dy 

Y 

= min J ( 32 ( 2 /, n(l + ipriTiyi + T 2 y 2 )))dy. 
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Let il’Tj be infinitely smooth functions approximating ipT in L^(0,Z). We will make a more 
precise choice of these fnnctions later. In order to incorporate the corrector ipTj in the equation of 
the approximating surface we need to understand the relation between the coordinates (t, s), x' and 
the underlying erF-periodic lattice. From x' = T{t)+sN = r{0)+tT+sN we derive t = T-{x'—T{0)). 
Applying the unfolding operator to ■ {x' — r(0))) we get 

= ilJTjiT ■ [e~^x' - e“^r(0)] + T -y) 

with y = {e~^x' — £:“^r(0)}, where [•] and {•} denote the integer and the fractional parts. Hence in 
order to align the argument with the lattice we need to shift it by adding t* = T ■ {e“^r(0)}: 

+ tl)) = i,TAT ■ - [e-AiO)]) + T-y)= ^tAT ' v) 

with y = {e~^x'}. On the domain [0, 1] of Fr we define a corrected tangent vector 

Te{t) := r(t) Y^l - (7.72) 

so that |Te(t)| = 1. For small values of e, the square root in the above is well defined. The 
approximating leading curve 7 ^ is given by integration of Tg, namely 

t 

le{t) ■■= j Te, 7 e ( 0 ) = 7 ( 0 ). 

0 

Now we set 

Ue(r(t) + sN) := 7 e(t) + sv, whenever r(t) + sA^ € 

Let us estimate the difference between u and Ug. First, it is clear that 

Te = T + 0(e), ne = n + 0(e), 7£=7 + 0(e), as e —>■ 0, (7.73) 

where := A v. Here and until the end of this section o(e), 0(e) are understood in terms of the 
L°“-norm. Hence 

llrte — l|L°o((.d[rr]) Oj as e —> 0 . 

Second, we have (c/. (j6.52p l 

V'us = (TiTe - T2 v\ T2Ts + Tiv) . 

It follows immediately that 

(7.74) 

We carry out this procedure for all rational arms ccPF^] leaving the surface unchanged on the rest 
of oj. We then assemble a new surface denoting it by in the same way it was done in the proof of 
Theorem IG. 81 using rigid motions. It is clear from the construction that 

Ug in in H^A)- 

Now we “add thickness” to the plate and define the approximating deformation as 

x^ 

usj^h'■= Usix') + hxsTisix') + A-^A(x',s ^x'), (7.76) 

where dj(x',y) is some function from C°°(uJ x Y) periodic with respect to y, which will be chosen 
later and the normal rig is defined by rig = diu^ A d 2 Ue on non-rational arms and bodies. The 
gradient of uSj,h is given by 

2;2 

^hUg-^h = (V'uel Ue) + hxA'^'us] dj) + A-^{V'dj + | 0). 
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Clearly, {Vu^, n^) € S0(3), thus we have 

h~^({V'ue\ne)^yhus^,h = Xsiy'uel Ue)^{Vnel dj) + o{e). 

On each arm a;(r) which does not belong to the set of rational arms u;(rj.) the surface itg is a rigid 
motion of u. Such a transformation preserves the matrix of the second fundamental form, i.e. for 
the matrix of the second fundamental form of we have Ilg = (V'ue)^V'ne = II on a;(r). 

Let us consider a rational arm a;(rr). Notice that for which is given via r' = Ken-si after a 
simple calculation using (17.721) and (I7.73P one has 


Ke = K + ^t)+0{e). 


Since n'^ = —KgTe we have 

diUe = —KeTeTi, Z = I, 2. 
Then from (I7.73p . (17.7411 and (17.7711 we obtain 


{VueyVne = -{k + + 0(e). 


(7.77) 


We conclude that 


h ^((V'ue, -I 3 ) 


a; 3 (ll| dj) + o(e) 

X 3 {- {k + K'ip'rl, j)ll\ dj) + 0{e) 


if x' € uiir n 5^, 
if x' € UJri 


(7.78) 


where dj := (V'rt, n)^dj, and uj = uirU U!ir U Ss is the natural decomposition of to into the union of 
rational arms ojr = (Jr the residual set of rational segments Ss and the rest of the set oJir- 

We understand the notation ijjTj as a family of functions whose members are related to rational 
arms via the corresponding directions T. 

Now we are ready to define a recovery sequence and prove the convergence of the elastic energy. 
We use the frame indifference of W, its Taylor expansion near I 3 , the unfolding operator and (17.781) 
to get 


h '^Eh{usj,h) = h '^Jw{e ^x',Vus^^h)dx = h ^ J W{e ^x',{y'ue\neyVhUSj,h)dx = 

= ^ y J Q 3 {y,%(iMdj))dx'dy + ^ J J Q 3 {y,Tei-iK +K'ipT^^)ll\dj))dx’dy + j 0{e)dx 

U}irC\S5 Y (jJr Y n 

+ ^ J Q3{£~^x',h~^{{\7'ue, ne)^yhUs,j,h-d:3))dx'. 

Notice that we can use the Taylor expansion of W due to the fact that II G L°°{uj) by Theorem 
16.81 Passing to the limit in the above as h —>■ 0 and taking into account the properties of two-scale 
convergence we arrive at 

^^^h~'^Eh{us^,h) = ^ J J Q 3 {y,ill\dj))dx'dy 

UlirClSs Y 

+ ^ J J Q 3 {y,{-{K +■ y))U\dj))dx'dy. 

LiJ-p Y 
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Choosing ipTj and dj, i.e. close enough to the minimisers of the problems (17.7011 and (14.221) with 
appropriate arguments we can make the term 

^ J J Q 3 {y,{-{K + ■ y))ll\dj))dx'dy 

OJir Y CJr Y 

arbitrary close to 

^ J QZmirn^'- 

(jJirUoJr 

For the term corresponding to the residual set Ss^ we have, in general, 

J Q 3 {y,{ll\d,))dx'dy> J QZ^ill)dx', 

Ssj xY 

however thanks to the estimate (17.711) it can be controlled. Thus, for each j we can choose ipT^j and 
dj so that 

\jim h-^Eh{usj,h) - Kom(w5j)| < -■ 

'h —>0 J 

The latter, together with (I7.71L implies that we can choose a subsequence hj —>■ 0 so that 

\h~^Eh{us^,hj) - ^^hom(w)| < -j- 

Passing to the limit as j —>■ oo we get the convergence of energies. Thus 

UT ■= USj,hj 

is a recovery sequence. The convergence of the gradients follows from (16.551) (17.751) and the con¬ 
struction procedure. □ 


A Appendix: Two-scale convergence and its properties 


Let w C be a bounded domain and Y := [0,1)"^ (in the present paper d = 2). We denote by [x] 
the integer part of x, and by {a;} its fractional part, i.e. [x] G {a:} := x— [a;] G Y. The definition 
uses the unfolding operator which maps functions of one variable x to functions of two variables x 
and y, where the “fast” variable y is responsible for the function behaviour on e-scale: 


Te{v){x,y) 


V (e [|] -I- ey) if a; G e(^ -|- T) C a; for ^ G 
0 otherwise, 


for a measurable on w function v(x). Notice that %(v)(x, y) is set to be zero on the cells intersecting 
with the boundary of w. This implies, in particular, that if n G then Te{v) G H^{Y)). 

One has 


vdx = 


Te{v)dydx Y J vdx, 


tjJ ujxY 


Ae 


where := {x G e(^ -I- T) O a; : e(^ + Y)r\duj ^ 0, ^ G Z'^}. 


Definition A.l. We say that a bounded in sequence converges weakly (strongly) two-scale 

2 2 

to V € LP{uj X Y) as e 0 and denote this by v^ ^ v (ve ^ v) if Te{ve) converges to v weakly 
(strongly) in LP{u! x Y). 


32 














We refer the reader to [5] for a comprehensive analysis of the properties of the unfolding operator 
7 ^, which, in particular, links the convergence of sequences in to convergence of their unfoldings 

in LF'{uj X Y). 

We also need a formal definition of two-scale convergence when not all variables have their fast 
scale counterparts, which is just a variation of the standard definition. 

Definition A.2. Using the notation accepted in this paper (i.e. il = ui x I, to C Y = [0,1)^, 
X = {x'jXs)), we say that a bounded in LP(D,) sequence Vg converges weakly (strongly) two-scale to 
V G X y) as e —?► 0 and denote this by ^ v (vg ^ v) if Te{ve) (where Te acts only with 

respect to x') converges to v weakly (strongly) in L^iQ. x Y). 

B Appendix: Formal asymptotic expansion for a recovery 
sequence 

In this appendix we present a formal asymptotic expansion for a recovery sequence which leads to 
the zero-determinant condition (15.4011 . 

Let u{x') be an arbitrary isometric surface, u € and -iplx^y) be a sufficiently smooth 

scalar function. As a starting point we consider the surface 

:= u{x') - e^-ipix', e-^x')n{x') (2.79) 

(similar to the case e ^ h^ e^). A simple calculation (c/. (14.3711 1 shows that (V'{tgP’^)^V'u|P’^ = 
2^2 + O(e^) at least formally. Ideally we need to add a sufficient number of corrector terms to 
the expression in (I2.79|) to make the error to be of order o{h). Let us consider a third order 
approximation, 

it^P’3 := u{x') — e^-iplx', e~^x')n{x') -\- £^ip{x', £~^x'), 
where is a sufficiently smooth vector-valued function. Taking the gradient we get 
V'u|P’3 = V'u — eVy-ijj n — (VV n + V^V'n -|- Vyip) 

Noticing that 

din = IIiii?i 4- IIi22?2i d2n = IIi27?i -I- Il222?2 

we have 

|aiM^P’3|2 = I _ 2^11^^ + 2Ri ■ dy^ip) Y 0(e3), 

|a2<P’3|2 = l+£‘^{{dy^ijjf - 21(1122 + 2i?2 • dy^pf) + 0(e3), 
diul^'^ ■ d2ul^’^ = £^[dy^-ipdy^-ip - 2iplli2 + Ri ■ dy^tpY R 2 ■ dy^pP) + 0(e3). 

Let us equate the second-order terms to zero. 

{dy,Pf - 2pilll + 2i?l ■dy,P> = 0, 

{dy^lpY - 2lpll22 + 2 .R 2 ■ dy^P> = 0, 
dyipdy^ip - 2^IIi2 -I- Ri ■ dy^p} + R 2 ■ dy^pj = 0. 

Clearly, the function ip has to satisfy certain solvability condition. Indeed, differentiating the above 
equations twice with respect to y, we get 

-Rl • 9yiV2V2^ = -i^mV2pf - 

^2 ■ dl,y,y,Y = -{dl,y,Pf " 9y,P + dl^y^pil22, 

dl,V^Pdv2V2p + id^y^pf + dy,Pdl^y^y^lP + dy^P dl^y^y^P - 2dl^y^ipili2 

YRi • dy^y^y^p} ^ 2 ’ ^y^y-^y^P^ — 0. 
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Eliminating the terms that contain derivatives of Lp and recalling that det II = 0 yields 


det(II + = 0. 


Acknowledgements 

This work was carried out under the financial support of the Engineering and Physical Sciences 
Research Council (Grants EP/H028587/1 “Rigorous derivation of moderate and high-contrast non¬ 
linear composite plate theories” and EP/L018802/1 “Mathematical foundations of metamaterials: 
homogenisation, dissipation and operator theory”). We would like to thank Dr Shane Cooper for 
helpful remarks. 


References 

[1] G. Allaire, 1992. Homogenisation and two-scale convergence, SIAM J. Math. Analysis, 23(6) 
1482-1518. 

[2] M. Cherdantsev, K. D. Cherednichenko, 2012. Two-scale P-convergence of integral functionals 
and its application to homogenisation of nonlinear high-contrast periodic composites, Archive 
for Rational Mechanics and Analysis, 204 445-478. 

[3] Cherdantsev, M., Cherednichenko, K.D., Neukamm, S., 2013. Homogenisation in hnite elas¬ 
ticity for composites with a high contrast in the vicinity of rigid-body motions. Preprint 
arXiv:1303.1224 

[4] P. G. Ciarlet, 1988. Mathematical Elasticity, Vol. I: Three-Dimensional Elasticity, North- 
Holland. 

[5] D. Cioranescu, A. Damlamian, G. Griso, 2008. The periodic unfolding method in homogenization. 
SIAM J. Math. Analysis, 40(4) 1585-1620. 

[6] E. De Giorgi, T. Franzoni, 1975. Su un tipo di convergenza variazionale. Atti Accad. Naz. Lined 
Rend. Cl. Sci. Mat. Natur., 58 842-850. 

[7] B. Dacorogna, 2000. Direct Methods in Calculus of Variations, Springer. 

[8] G. Friesecke, R. D. James, S. Muller, 2002. A theorem on geometric rigidity and the derivation 
of nonlinear plate theory from three-dimensional Elasticity. Comm. Pure Appl. Math., 55(11) 
1461-1506. 

[9] G. Friesecke, R. D. James, S. Muller, 2006. A hierarchy of plate models derived from nonlinear 
elasticity by Gamma-convergence. Arch. Rat. Mech. Analysis, 180(2) 183-236. 

[10] G. Griso, 2004. Error estimate and unfolding for periodic homogenization. Asymptot. Anal., 
40(3-4) 269-286. 

[11] G. H. Hardy, J. E. Littlewood, G. Polya, 1952. Inequalities. Cambridge Mathematical Library 
(2nd ed.), Cambridge University Press. 

[12] P. Hornung, 2008. Approximating isometric immersions. Comptes Rendus Mathematique, 
346(3) 189-192. 

[13] P. Hornung, 2011. Approximation of flat isometric immersions by smooth ones. Archive 
for Rational Mechanics and Analysis, 199(3) 1015-1067. 


34 


[14] P. Hornung, S. Neukamm, I. Velde, 2014. Derivation of a homogenized nonlinear plate theory 
from 3d elasticity. Calculus of Variations and Partial Differential Equations, 51(3-4) 677-699. 

[15] S. Muller, M. R. Pakzad, Regularity properties of isometric immersions, 2005. Mathematische 
Zeitschrift, 251(2) 313-331. 

[16] S. Neukamm, 2010. Homogenization, linearization and dimension reduction in elasticity with 
variational methods, Ph.D. Thesis, Munich. 

[17] S. Neukamm, H. Olbermann, 2014. Homogenization of the nonlinear bending theory for plates. 
Calculus of Variations and Partial Diff. Eq.. 

[18] G. Nguetseng, 1989. A general convergence result for a functional related to the theory of 
homogenization, SIAM J. Math. Anal, 20 608-629. 

[19] M. R. Pakzad, 2004. On the Sobolev space of isometric immersions, J. Differential Ceom., 66(1) 
47-69. 

[20] O. Pantz, 2001. Quelques problemes de modelisation en elasticite non lineaire (Some modeling 
problems in nonlinear elasticity). These de I’Universite Pierre et Marie Curie. 

[21] M. Spivak, 1979. A Comprehensive Introduction to Differential Geometry, Vol. Ill, Publish or 
Perish Inc. 

[22] I. Velcic, 2013. A note on the derivation of homogenized bending plate model. Preprint 
arXiv:1212.2594. 

[23] V. V. Zikov, S. M. Kozlov, O. A. Oleinik, 1994. Homogenisation of Differential Operators and 
Integral Functionals, Springer. 

[24] V. V. Zhikov, S. E. Pastukhova, 2005. On operator estimates for some problems in homogeni¬ 
sation theory. Russ. J. Math. Phys., 12(4) 515-524. 


35 


